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A GENERAL SMALL-DEFLECTION THEORY FOR FLAT SANDWICH PLATES 

■ By Charles Libove and S. B. Batdobf 


SUMMARY 

^4 small-dejledion theory is developed Jor the elastic behavior 
of orfhotropie flat plates in which deflections due to shear are 
taken into account. In this theory, which covers all types of 
flat sandwich construction, a plate is characterized by seven 
physical constants (five stiffnesses and two Poisson ratios) of 
which six are independent. Both the energy expression and 
the differential stations are developed. Boundary conditions 
corresponding to simply supported, clamped, and elastically 
restrained edges are cortsidered. 

INTRODUCTION 

The advent of high-speed flight and the concurrent neces- 
sity of maintaining aerodynamically smooth surfaces under 
high stress have led to the increased study of sandvdch-plate 
construction as a possible substitute for sheet-stringer con- 
struction in airplane design. A sandwich plate consists 
essentially of a relatively thick, low-density, low-stiffness 
core bonded between two thin sheets of high-stiffhess ma- 
teriaL Materials that have been considered for the core 
include balsa wood, hard foam rubber, cellulose acetate, 
resin-impregnated cloth fashioned into a honeycomb, 
corrugated metal sheet, and even closely spaced- stiffeners 
of the conventional typo. The face sheets may be of metal, 
ply^vood, wood-pulp plastic, or some other type of high- 
stiSness material. 

Because of the low-stiffness core, the sand-wich plate will, 
in general, experience appreciable deflection due to shear. 
Furthermore, becaiise the face sheets or core (or both) may 
have orthotropic stretching properties, the sandwich plate 
will in general be orthotropic in its flexural properties. As 
a result, ordinary plate theory, which is based on the assump- 
tions that the plate is isotropic and that deflections due to 
shear may be neglected, cannot be used to determine the 
stresses, deflections, or buckling loads of sandwich plates. 

A general small-deflection theor}’- for flat orthotropic plates 
is therefore developed in which deflections due to shear are 
taken into accovmt. The theory is applicable to any type 
of orthotropic or isotropic sandwich that behaves essentially 
as a plate, pro-vided certain physical constants are kno-wn. 
These physical constants (two flexural stiffnesses, two shear 
stiffnesses, a t-wisting stiffness, and two Poisson ratios defined 
in terms of curvatures) serve to describe the plate deforma- 
tions associated Avith simple loading conditions and may be 
regarded as fundamental properties of the plate. For 
simpler types of sand-wich construction the physical constants 
can be evaluated theoretically from the geometry and physi- 
cal properties of the materials used. For more complicated 
types of construction, these constants can be evaluated by 
means of simple tests on samples of the assembled sand-wich. 


as described in appendix A. A reciprocal relationship 
betw'een the flexural stiffnesses and Poisson ratios is derived 
in appendix B. 

As is the case -with ordinary plate theory, the orthotropic 
plate theory consists of two parts, each complete in itself. 
These parts are a set of six differential equations, three of 
which express the equilibrium of an infinitesimal plate ele- 
ment and three of which relate the curv^atures and twist of 
the clement to the forces and moments acting upon it, and 
an expression for the total potential energy of the system 
comprising the plate and the forces acting upon it. The six 
differential equations involve six variables. However, it is 
sho-wn how these simultaneous equations can be reduced to 
a single equation of sixth order involving any one of the 
variables alone. In appendix C the consistency between the 
differential equations and the potential-energy expression is 
sho-wn by a variational method. 

The consideration of deflections due to shear makes neces- 
sary the specification of one more boundary condition than 
in ordinary plate theory. T^iis fact was first appreciated by 
Reissner in reference 1. Because of some arbitrariness in the 
choice of the additional boundary condition, two types of 
simple support and two types of clamped edges are possible. 
Furthermore, three boundary conditions can be specified 
for a free edge, in contrast to ordinary plate theory. Bound- 
ary conditions more general than freedom, simple support, 
or clamping are considered in appendix C. 

A number of investigations related to the problem of ortho- 
tropic- or isotropic-sand-wich-plate analysis have been made 
pre-viously. Theories for the bending of orthotropic plates 
due to lateral loads and budding due to edge loads, neglecting 
deflections due to shear, are given in references 2, 3, and 4 
and pages 380-384 of reference 5. The effect of shear on 
the bending due to lateral load of homogeneous isotropic 
plates and isotropic sand-nich plates is considered in refer- 
ence 6. The effect of shear on the bending due to uniform 
lateral load and buckling due to edge compression of simply 
supported isotropic sandwich plates with homogeneous cores 
is considered in investigations by Hopkins and Pearson and 
by Leggett and Hopkins. A rough method of taking into 
account deflections due to shear in the buckling of simply 
supported orthotropic sandwich plates is used in reference 7. 

The present theory may be regarded as a natural exten- 
sion to plates of the approximate theory used in pages 170- 
174 of reference 8 to take into account deflections due to 
shear in a beam. The theory of this paper is more general 
than the aforementioned theories in that it applies to ortho- 
tropic or isotropic sandwich plates -with homogeneous or 
nonhomogeneous cores and -with arbitrary boundary condi- 
tions, it presents both the differential equations and the 
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energy expression for the plate, and it is applicable to prob- 
lems that involve lateral as well as edge loads.The differential 
equations of the present theory are reduced to special forms 
in order tliat they may he 'compared with the . equations 
obtained in references 6 and 6. 

The detailed development of the theory comprises most 
of the following sections and the appendixes. The main 
parts of the theory are summarized briefly in a section entitled 
“Recapitulation of Principal Results.” 


SYMBOLS 

z, y, z orthogonal coordinates; z measured normal to 
plane of plate and z an(L2/ parallel to principal 
axes of flexural symmetry, inches 
w deflection of middle surface of plate, measured 

in 2-direction, inches _ _ . 

q intensity of lateral loading, pounds per square 

inch 

mtensity of internal shear_acting in 2-direction 
in a cross section originally parallelto.t/2-plane, 
pounds per inch 

Qg intensity of internal shear acting in 2-direction in 

a cross section originally parallel to rz-plane, 
pounds per inch 

Mx intensity of internal bending moment acting 

upon a cross section mriginaUy parallel to 
3/2-plane, inch-pounds per inch 
mtensity of internal bending moment acting 
upon a cross section originally parallel to 
ss-plane, inch-poimds per inch 
Mjy intensity of internal twisting moment acting in a 

cross section ox'iginally paiallel to y2-plane or 
ar2-plane, inch-pounds pOT inch 


N, 


N, 


Dx, 


-O — 

I^Bendin^ 


Nx intensity of middle-plane tensUe force parallel to 

r2-plane, pounds per inch 
intensity of middle-plane tensile force parallel to 
y2-plane, pounds per inch 
intensity of middle-plane shearing force parallel 
to yz-plane and rz-plane, pounds per inch 
flexural stiffnesses of plate with anticlastic 
bending unrestrained, inch-pounds 
moment per incIA 
Curvature / 

twisting stiffness of plate, inch-pounds 
/ Twisting moment per inch '\ 

\ Twist / 

flexural stiffness of ordinary plate, inch-pounds 
shear stiffnesses of plate, pounds per inch 
Poisson ratios for plate, defined in terms of 
curvatures 

Poisson ratio for ordinary plate 
shear-strain angles due to shears Qx and Qn, 
respectively, radians 
h thickness of plate, mehes 

a, b length and width, respectively, of rectangular 

plate, inches 

V total po tential energy of system, mch-pounds 

' strain energy of bending of plate, inch-pounds 
T^a - potential energy of external loads, inch-pounds 
u, V displacements in a:-direction and y-dircction, 

respectively, of a point in middle surface of 
plate, inches 

[Z7], M, [A^, [P] differential operatora 

SIGN CONVENTION 


D: 


D - 

Dq^, Dq^ 
fix, /ty 


yx,yt 


The sign convention and notation used in the present 
paper are, wherever convenient, tlie same as those used by 
Timoshenko in reference 5. 



Fiqube 1.— Forces and moments acting on dldeientlal element dx dt. 
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The X-, y-, and 2 -axes of an orthogonal coordinate system 
are oriented so that the zy-plane coincides \rith the undis- 
torted middle plane of the plate. Deflections w are measured 
normal to the ary-plane and are positive in the positive direc- 
tion of the 2 -axis. The lateral load q is also positive m the 
direction of the 2 -axis. 

The internal shears and Qy, moments M^, and M*),, 
and middle-plane forces Nz, N,, and iVi, are shown m figui-e 1 
acting in their positive directions upon an infinitesimal 
element of length dx and width dy cut from the imloaded 
plate by planes parallel to the a: 2 - and ye-planes. Only the 
forces and moments acting^ on two adjacent faces of the 
element are shown. The forces and moments on the oppo- 
site faces differ from those on the faces shown only by infini- 
tesimal amounts. The directions in which they act, however, 
are opposite (for example, moment A/* dy on the face shown 
is counterclockwise; moment A/* dy on the opposite face 
would be shown acting clockwise). The twisting moment 
and middle-plane shearing force acting on any cross section 
are known, from equilibrium considerations, to be equal to 
the twisting moment and middle-plane shearing force acting 
on a cross section at right angles. The symbols Mxy and 
therefore appear in both of the faces shown in figure 1 . 

For convenience, in this report the 2 -direction is sometimes 
referred to as the vertical direction and planes parallel to 
the a^-plane are sometimes referred to as horizontal planes. 

PHYSICAL CONSTANTS 


The physical properties of the plate are described by means 
of seven constants: the fle.xural stiffnesses D* and the tvdst- 
ing stiffness the transverse shear stiffnesses and 
and the Poisson ratios Hx and ju„. Definitions of these 
constants are obtained by considering the distortions of the 
differential element of figure 1 under simple loading conditions. 

Let all forces and moments acting on the element be zero, 
except for the moments Mx acting on two opposite faces. The 

effect of Mx is to produce a primary curvature in the 

middle surface of the element and also a secondary curvature 

which is a Poisson effect. Then Dx is defined as the 

negative of the ratio of moment to primary curvature or 


D.= - 


Sa:* 


( 1 ) 


when only Mx is acting, and fix is defined as the negative of 
the ratio of Poisson curvature to primary curvature or 


u ^ 


( 2 ) 


when only Mx is acting. No other distortions are assumed 
but and when Mx acts. The minus signs are intro- 
duced in order to make Dx and nx essentially positive 
quantities. 


Similarly, D, and n, are defined as 


I/O 

1 

II 

Cl 

C3) 



b^V) 


bi? 

(4) 

&2/* 



when only A/, is acting. 


If, now, aU of the forces and moments are equal to zero 
except Mx, acting on aU four, faces, the only distortion 

produced is a twist and Dx, is defined as the ratio of 

twisting moment to twist or 


D 


bxby 


(5) 


when only Mx, is acting. 

The transverse shear stiffness Dq^ is defined by letting only 


the shears Qx act on opposite faces of the element (except for 
an infinitesimal moment of magnitude Qx dy dx required for 
equfiibrium). The distortion is assumed for the moment 
to be essentially a sliding of one face of the element with 
respect to the opposite face, both faces remaining plane. As 
a result of this shdmg, the two faces parallel to the a; 2 -plane 
are distorted from their rectangular shape into parallelograms 
by an amount y*, which is the shear angle measm-ed in the 
arz-plane. The shear stiffness Dq^ is defined as the ratio of 
shear to shear angle or 

( 6 ) 


when only Qx is acting. If the sides of the element are kept 
parallel to the 2 -axis, the slope of the middle surface is 


bx 2^ 


when only Qx is acting. 

In a similar manner, the shear stiffness Dq^ is defined as 
the ratio of the shear on the faces parallel to the s 2 -plane to 
the shear angle measured in the y 2 -plane when only Q, is 
acting or 

m 


when only Q, is acting. If aU sides of the element are kept 
parallel to the 2 -axis, the slope produced is 




when only Q, is acting. 


The constants just discussed serve to define the ortho tropic 
sandwich plate; they can be evaluated theoretically if the 
properties of the component parts of the sandwich are known 
and if the plate is of simple construction. In any event, the 
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coustants can be determined experimentally by means of 
bending tests and twisting tests on beams and panels of the 
same sandwich construction as the plate. A description of 
the tests required is given in appendix A. 

Although seven physical constants. have been, discussed, 
they need not.all be independently determined for if any three 
of the four constants Z>*, Dy, ft*, and y.y are known the fourth 
can be evaluated from tlie relationshp 

HrDy = [iyDx _ ( 8 ) 

This relationship, based on a generalization of iSlaxwell’s. 
reciprocal law, is derived in appendix B. 

The shear stiffnesses Dq^ and merit some additional 
discussion. The distortion due to shear was assumed to be 
a sliding of the cross, sections over each other, the cross secK 
tions remaining plane and the shear strains remaining con- 
stant for the entire thickness of the plate and equal to the 
shear angle y* or y,. Actually, if the plate is continuous 
enough for cross sections to exist at all, under shear the cross 
sections generally tend to warp out of their plane condition 
(p. 170 of reference 8); ths warping makes the shear angle, 
as defined for equations (6) and (7), meaningless. The shear 
strain varies with depth and an average shear strain will have 
to be used as the effective shear angle y* or y„ for purposes 
of defining effective shear stiffness or Dq^. If the exper- 
imental method is used . (see- appendix A), this difficulty is 
not encountered because, instead of a shear angle, curvatures 
are measured, and the stiffnesses, obtained are automatically 
the effective stiffnesses. 

Despite the general tendency of cross sections under shear 
to warp, the assumption that they remain plane (though not 
normal to the middle surface) can be. shown .to. be almost, 
correct for those sandwiches in which the stiffness of the 
core is very small compared with tlie stiffness of the faces 
(for example, Metalite, honeycomb). For such sandwiches 
the shear stiffnesses Dq^ and Dq^ can. be readily calculated, 
because the faces may be assumed to take all the direct bend- 
ing stress and the vertical shear may therefore be assumed 
uniformly distributed in the core. The shear angles y* and 
y„ will then be constant thi-oughout the core. 

For those sandwiches in which cross sections under shear 
may not he assumed to remain plane, the tendency of these 
cross sections to warp introduces a further, complication 
which can, however, be resolved by means of a justifiable 
simplifying assumption. This complication is due to the 
fact that if the cross-sectional warping is partially or com- 
pletely prevented the effect will be to increase, the shear 
stiffness Dq^ or Dq^. The shear stiffnesses, thus, depend not 
only on the properties of the plate materials but also on the. 
degree of restraint against cross-sectional warping. For the 
purpose of the present theory the shear stiffnesses Dq^ and 
Dq^ are assumed to be constant throughout the plate and 
have the values they would have if cross sections were allowed 
to warp freely. The error caused by this assumption will 
be mainly local in character, being most pronounced in the 
region of a concentrated lateral load, where a sudden change 
in the shear tends to produce a sudden change in the degree 


of. warping which is prevented by continuity of the plate. 
The error will probably be negligible in the case of distributed 
loads, for which there are only gradual changes hi the shear. 
A discussion of this error, in connection with beams is con- 
tained in pages 173-174 ot.referenee. 8 and.in reference 9, 


DIFFERENTIAL EQUATIONS FOR PLATE 
DISTORTION EQUATIONS 


5’w 


Equations can be derived relating the curvatures and 

and the twist fl't any point in the plate to the inter- 
nal shears and moments acting at that point. 

Equation for the curvature — An expression can be 

obtained for the total cm-vature hi the x-dircction by 

adding together the contributions made by each of the shears 
and moments acting separately. From equation (1) the 
curvature contributed by Mx is found to bo 

Mx 

~lTx 

Equations (3) and (4) can be solved for the contribution to. 
by My which is 

D, 

Finally, the equation following equation (6) indicates that 
the existance of produces a curvature in the middle 
plane eq^ual to 


1 ^ 
^ i>x 


The moment M*„ and the shear Qy make no contribution to 

* 

Addition of the three component curvatures gives 
Mx 1 My I 1 iiQx 


Equation for the curvature — Similar considerations 


(9b) 


give the ciu-vature in the i/-direction as 

My. Mx, 1 ^Qv 

Equation for the twist expression for the twist 

is obtained by firat writhig an expression for the twist- 
ing moment Mxy hi terms of the distortions of the clement 
dx dy. 

Let the middle surface of the element bo distorted so that 

it acquires a twist Furthei- assume that each line 

elementnoi-mal to this middle surface before its distortion (a) 
fh-st rotates .so as to remain normal to it afttT its distortion, 
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(b) tben rotates through an angle Yi in a plane parallel to 
the as-plane, and (c) then rotates through an angle y, in a 
plane parallel to the i/a-plane. (Rotations (b) and (c) pl*o- 
duce parallelogram-type distortions of cross sections and are 
therefore denoted as shear angles Yi and y„.) 

Distortion of the element as a result of rotations of type 
(a) is shown in figure 2 (a). Distortion of the element due 
to rotations of type (b) is shown in figure 2 (b) on the assump- 
tion that Yi is zero at the center of the element and is chang- 
ing imiformly in the y-direction. Distortion of the element 
due to rotation of type (c) is shown in figure 2 (c) on the 
assumption that is zero at the center of the element and 
changing uniformly in the a:-direction. The magnitudes of 
the displacements shown in figure 2 are obtained by consider- 
ations of geometry, the details of which are not given. 

The t%visting moment acting on aU four cross sections 
of the differential element is proportional to the, shear strain 
of the upper and lower surfaces, because this type of strain, 
throughout the thickness of the element, produces the hori- 
zontal shearing couples that make up By superposi- 

tion of the three distortions shown in figure 2, the shear 
strain in the upper (or lower) surface can be written as 


2 bxby 


7t h 5yi ^ &Ytf 

2 iixbyj 2 dy 2 dr 


/ yw 1 dYi 1 &yA 

\bxdy 2 by 2 dx J 


5/ 



sorrace. 

Fiouek 2. — Distortions of element t£r (fr in twisting. 


and theretore, 


AlxgOZh 


bxby 


1 

2 by 


2 bx J 


Substitution for yi and y» id terms of Qx and Qf (equations 
(6) and (7)) gives . 


Mx,=h' 


'd^w 1 1 dQx 1 1 bQ,\ 
bxby 2 dy 2 Dq^ dr J 


where A' is a proportionality constant absorbing A. When 
Qx and Qg are both set equal to zero, the above equation 
must reduce to equation (5), because only il/,, is acting on 
the differential element. The constant A' is therefore iden- 
tified as D„ and the equation for twisting moment becomes 



drdy 


i 

2Dq^ by 


1 J_^\ 

2 Dq^ dr J 


Solution for 


b‘w 

drdy 


yields the following equation analogous to 


d% 


the equations already obtained for and 


d®ia. 

W' 


b‘w Alxu ,1 1 'dQx.l 1 ^ 
drdy Dq^ by "^2 dr 


(9c) 


EQUUilBBIUM EQUATIONS 

The element dx dy must be in equilibrium under all the 
forces and moments acting upon it. This condition implies 
that certain relationships must exist among these forces and 
moments. These relationships can be derived by considering 
the changes that occur in the forces and moments from one 
face to the opposite and writing the equations of equilibrium 
for the element. The equations are the same as in ordinary 
plate theory. For equilibrium of forces in the r- and 
y-directions, these equations are obtained from equations 
(196) of reference 5 : 


dzV, bNxx „ 
dr ‘ dy 

(10a) 

o 

II 

ClOb) 


The equation for equilibrium of vertical forces is given at the 
top of page 305 of reference 5 as 


b^Mx 


d*3/, - / 

dr* 

drdy 

^ by* V 


{i+Nx 


(11a) 
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And the equations for equilibrium of moments about the y- 
and a^-axes are obtained from equations (188) and (189) of 
reference 5 as 


dy dr 


(11b) 




dr 


dMy 

dy 


(11c) 


(Equations (11) are also derived in appendix C by minimiza- 
tion of the potential energy.) Note that the left-hand side 
of equation (11a) can, by virtue of equations (11b) and (11c), 
be simplified to 

I ^Qv 

dx dy 

If, as is customary in small-deflection theory, the middle-, 
plane stresses iV*, Ny, and Nxy are a^umed to be unchanged 
in the course of the plate’s deflection and equal to their initial 
values before application of lateral load, then equations (10) 
are automatically satisfied and equations (9) and (11) con- 
stitute the six fundamental differential equations that deter- 
mine the .forces, moments, and distortions throughout, the 


orthotropic plate. They can be used in their present form or 
in the alternate-form obtained in the following section. 

ALTERNATE FORM OF THE DIFFERENTIAL EQUATIONS 

The fundamental differential equations (9) and (11) can be 
transformed so as to separate variables. Equations (9) are 
first solved for M*, My, and Mxy to obtain 



With the left-hand side of equation (11a) simplified to 
and the above expressions for A/*, My, and 

substituted into equations (lib) and (11c), equations (11) 
become, after some regrouping of tei’ins, 


5S^)»+(|) «.+(|) e.=-s 


r r. , ^x] , riDxy y , Dx b' ,1^ ,n 

£)» 


~| n _0 

Do bxdy^ {l-iXxiii)Dn bxdijJ 


r_z? ^ 

L ^^•'dx^dy I-Mx/xA' 


d» I y 
dx‘dy~^dy’ 








dxdy 


X o , fl Px, y I D. S' _ t) -n 

,)Do, 5xdyJ L2 ^y dx^~>~ M 


These three equations can be solved to obtain a differential equation for w alone in terms of g, an equation for Q* alone in 
terms of q, and an equation for Qy alone in terms of g. This separation is accomplished most easily, for the case in which 
JV*, Ny, and Nxy are constant throughout the plate, by treating the tliree differential equations as though they were alge- 
braic equations and solving for iv, Qx, and Qy by means of determinants. The terms in the determinants arc the differential- 
operator coefficients of w, Qx, and Qy appearing in the three equations. In expanding these determinants, the rule for multi- 
plication of hnear operators must be used. For example. 


d’ d* 

dp bxdy dxdp 


As a result of such a solution, the following differential equations are obtained for w, Qx, and Qyi 


[D]w=—\M\q 


[D]Q.= -[N]q 
[D\Qy^-[F\q 


(13a) 

(13b) 

(13c) 
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where [Z?], \hi\, [A^, and [P] are differential operators defined as 


1 DxyDx b® . j 

^ J J-J DxDy 2 BxyDxHy 2 BxyDyfXx 

1 l( 

^ \ J) J) BJDy ^DxyDxfiy ^DxyDyll^ 


2 Dq^ dx®'' 

iv2 flo. + J 

bifby^ ’ 

1,2 + T}7T j 

' bx‘by* 


\ So^S^; / 2>x-*5W+2 15^ V 

ri DxrO-—t^t) I \-2N ■ ■^— ^— r~ f 

(1— (■^® 5rdy) 


+ 


baray' 


bx^iyy' 


7.+N,^+2N., 


a* 

a^aj/* 


’)+ 


(14a) 


= 215^^ 


^ DxDy cj DxyDxHy 2 DxyDyflx 

1 1 

1 DxyDy d* 

ri Z?*»(l— Mi/tn) 


\ DqDq^ J 

'bx^by^^ 

^2DQpQ^by^ 

l2 Dq^ 

^DoJ bp 


[^']=4 


[P]=4 


1D„D,V ®''°''‘'“5®'’®''‘''i 

1-2 &x»+\ z?^;; / 

^ DxDy ■g DxyDxliy 

\ K 


a» 


a® , 1 DxtDy S® 71 rr> \ I n 1 ^ 

a?a^+2 ba:2>2/« ^*5? 5^ 


1 D„Dg a* 


2 z?o aaj^ay 




~2 a® , 1 PxyPy &® 


a® 


a* 


ax®a2/>'^2 z?«/a2^ 


(14b) 

(14c) 

(14d) 


Equations (12) and (13) taken together constitute an alter- 
nate set of differential equations that the plate must satisfy. 

COMFABISONS WITH PEEVIOUS SOLUTIONS 

Homogeneous isotropic plates, deflections due to shear 
neglected. — The usual fourth-order equation for homogene- 
ous isotropic plates, in which deflections due to shear are 
neglected, can be obtained from equation (13a) by letting 

Dq=Dq= CO 


y^=y.^ = H 

Di=Dy=B{l—ii^ 
Diy=D0- — p) 


With these substitutions made, equation (13a) becomes, after 
some transposition of terms, 




5*10 5®tc_ 1 

b2*5y*'^&2/* ~D 


(<Z+JV, 


g+A^,0+2ATx 


IN} \ 
' bxdyj 


which is the same as equation (197) of reference 5. 

Isotropic sandwich plates, deflections due to shear 
considered. — ^The differential equations for isotropic sand^ch . 
plates are obtained in reference 6 by use of Castigliano’s 
theorem of least work for the case in which the middle- 
surface forces iV*, Ny, and iV„ are zero. The equilibrium 
differential equations of reference 6 are equivalent to equa- 
tions (11) of the present paper. Equations (10a), (lOd), 
(lOe), and (lOf) of reference 6 can be solved simultaneously 
to obtain the following equations for the curvatures and 
twist in terms of the vertical shears and moments (the 
notation is that of reference 6) : 

5*u> M. , j>My ,1 bVr 1 /5F, , bVy\ 

5r»- Z>(1-«-®)"^ZI(1 -f®)‘^( 7. dr VA ^ iyy J 

bho My , vMx . 1 1 fi>Vx . d^A 

dy® - Z?(l-i/®)"''Z?(l-j^ '^C, by C\bx^ by ) 

d*w , 1 (^Vx , bVy\ 

bxby Z?(l -p) ^ dy dx 
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The symbols H, r*, and 1"^ in the above equations correspond 
to — Jl/iy, Qx, and Q^, respectively, in the notation of the 
present paper. The quantities D, Gi, C„, and v are physical 
constants for the plate. The above equations are seen to be 
identical in form to equationsr (9). of the present paper (il Dq^ 
is set equal to for isotropy in the x- and y-directions) 

except for the additional cm^a- 


ture equation. This term arises from tlie consideration of 
stresses and strains in the vertical direction, which were neg- 
lected in the present paper on the ground that they have a 
negligible, effect on the over-all flexural behavior of the plate 
and ai‘6 only important in the neighborhood of concentrated 
loads. Setting Gb equal to infinity makes the equations 
derived from reference 6 completely identical in form to equa- 
tions (9) of the present paper. It should be mentioned that 
the quantity w as used in reference 6 is not the deflection of 
the middle surface but “a weighted average across the thick- 
ness of the deflections of all points of the plate which lie on a 
normal to the middle surface.” 


is applied only to the middle surface at the boundary and no 
horizontal forces are applied to prevent the y-displacement of 
other points in the boundary, then M,,, which is made up of 
such horizontal forces, must be zero. Two different types of 
simple support thus emerge. For simple support in which all 
points in the boundary are prevented from moving parallel 
to the edge, the conditions are 


w=0 

(ICa) 

o 

II 

(lOb) 

II 

o 

(IGc) 


For simple support in which all points in the boundary, 
except those in the middle surface, are free to move parallel 
to the edge, the conditions are 

10=0 (17a) 

M,=0 (17b) 


BOUNDARY CONDITIONS 

The boundary conditions are first discussed for those types 
of edge support most commonly assumed in practice namely, 
complete freedom, simple support, and clamping. (More 
general kinds of support are considered in appendix C.) 
These supports are characterized by the condition that no 
work is done by the moments and vertical forces at the 
boundary. A boundary parallel to the y-axis is considered ; 
the conditions for a boundary parallel to the r-axis can be 
obtained by replacing r by y and vice versa, except in the 
subscripts of Mxy and iV*„. 

Free edge. — ^The boundary conditions for a free unloaded 
edge parallel to the j/-axis express the conditions of zero 
bending moments, zero twisting moment, and zero vertical 
force, or 

Mx =0 .. . (15a) 

M„=0 . (15b) 

<2^=0 (15c) 


If the free edge carries load, the middle-plane forces iV, and 
Nxv vdU not in general be zero and the boundary condition 
of zero net vertical force becomes 




(15c') 


instead of equation (15c). 


Simply supported edge. — The principal boimdaiy con- 
ditions for a simply supported edge parallel to the y-axis are 
«!=0 and M*=0. If to these two conditions is added the 
restriction that there is . no ^-displacement of points in 
the boundary, then the shear angle yy is zero and therefore 
Q 

yy^=0. If, on the other hand, the support at the boundary 


Mxy=Q 


(17c) 


Of the two types of simple support, the first (equations 10)) 
is more: likely to occur in practice. 

Clamped edge. — The principal conditions characterizing 
a clamped edge parallel to the y-axis are zero deflection of the 
middle surface and zero rotation of the cross sections making 
up the boundary (that is, the boundary plane remains parallel 
to the z-a.xis). The requirement of zero deflection is satisfied 
by letting w=0 at the boundary. The requirement that 
boundary cross sections remain parallel to the z-axis is satis- 


fied by letting 


as the equation following equation (G) 


indicates. (Note that if deflections due to shear are 


neglected by letting Dq^= «>, then the last boundary condi- 


Qw 


tion reduces to ’§^=0) which is familiar in ordinary plate 
theory.) Just as in the case of simple support, the tim'd 
boundary condition is either -^=0 or Mty=Q depending on 


whether or not points in the boundary (other than those 
points in the middle surface) are prevented from moving 
parallel to the edge. Thus, two types of clamping are 
possible. For a clamped edge in wliich the points in the 
boundary of the plate are prevented from moving parallel 
to the edge, the conditions are 


v)=0 

(18a) 

II 

O 

(18b) 

Qv n 

(18c) 


For a clamped edge in which the points in the boundary 
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(except those in the middle surface) are free to move parallel 
to the edge, the conditions are 


w=0 

bx 

0 


(19a) 

(19b) 

(19c) 


The latter type of clamping is very unli kely to occur in prac- 
tice, because any practical type of restraint that keeps the 
boundary from rotating has to be applied over an appreciable 
part of the thickness of the edge and therefore prevents most 
points in the boundary from moving freely parallel to the 
edge. 

The boundary conditions just discussed, as well as bound- 
ary conditions corresponding to more general types of 
support, are derived in appendix C by a variational method. 

POTENTIAL-ENERGY EXPRESSION 

STBAIN ENEEGT 

An expression can be obtained for the strain energy T”! 
produced by the moments ilf*, and M„ and the shears 
Qg and Qy by considering the work done by these moments 
and shears in distortiug the differential element of figure 1. 

The work of the moments il/i dy is equal to ^ M^dy times 

the counterclockwise rotation of the right-hand face %vith 
respect to the left4iand face of the element. This rotation 
is made up of two parts; the rotation caused by the moment 
Mg itself and the Poisson rotation caused by the moment My. 

The sum of these two parts is — ^ dx. (Note 

that although the term makes a contribution to the 

cun'ature of the middle surface, this term represents a rate 
of change of sliding rather than a rate of change of rotation 
and therefore makes no contribution to the rotation of one 
face with respect to the opposite.) The work of the moment? 
Mr is therefore 

-i Mrdy ^ dx 


or 


I (Ml 
2\Dr 


MxMy \ J J 

Hy^^yixdy 


Similarly, the work of the momenta M„ is 


2 \ Dy 


MxMy 


^ dxdy 


( 20 ) 


( 21 ) 


The work of those moments Mxy acting in the faces parallel 
to the ars-plane is equal to ^ Mxy dx times, the clock%vise rota- 
tion of the nearer face (as seen in fig. 1) with respect to the 


farther face. This rotation is made up ot the two parts 
shown in figures 2(a) and 2(b) and is equal to 

bxby ^ &j/ ^ 

or, replacing y* by its equivalent in terms of Qx (equation (6)) , 

The work of the moments AJxr parallel to the arz-plane is 
therefore 






^ bQ, 

bxby Dq by 


dx dy 


Similarly, the work of those moments Mix parallel to the 
« 2 -plane is 

1 , , / b^w 1 7 >„ 

2^^^»{jbxby bx)^^"^^ 

The total work of the moments is, by adding the last 
two expressions, 


Mx 


(ihli 1 1 bQx 1 1 a<3A 7 7 

2 I5qj. by 2 Dq^ bx ) 

Mx 


The factor in parentheses is simply from the equation 

preceding equation (9c), and the work of the moments Mxy 
therefore becomes 

^dxdy (22) 


The work of the shears Qx '^s^Qxdy times the downward 

distance through which the right-hand face slides. _^th 
respect to the left-hand face. This distance is y* dx and 

work is therefore ^ Qxyx dx dy. Keplacement of yi by its 

equivalent in terms of Qx gives 




( 23 ) 


for the work of the shears Qx. Similarly, the work of the 
shears Qy is 

( 24 ) 

Integration of the energy expressions (20) to (24) oyer the 
entire plate gives, as the total strain energy due to bending 
and shear, 

’’■-5 1 1 






(25) 
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Elimination of M*, M^, and Mxy by use of equations (12) transforms the strain-energy expression (25) into 


F =i r r ^ FA <2x \ T , Q, XT 

^ 2j J (1— /x*/Hj,L£»a:\ba: A ^ 1— AtxMt- L^a:\£te l^/JL^yW 



In addition to the strain enei'gy of bending and shear, there is the energy of stretohing of the middle surface produced 
by the forces N^, Ny, and Nxy. In small-deflection theory, these forces are assumed to remain constant during lateral deflec- 
tion. The strain energy of middle-surface stretching is therefore a constant independent of the lateral deflection. This 
energy does not affect any solution and may be omitted from consideration. 

POTENTIAL ENERGY OF EXTERNAL FORCES 

The potential energy acquired by the external forces in the course of the lateral deflection of the plate is independent 
of the internal construction details of the plate and depends only on the displacements of the middle sui-face. The potential- 
energy expression for the orthotropic sandwich plate is therefore the same as for the ordinary homogeneous isotropic plate; that 
part of the expression due to the forces Nx, Ny, and Nxy at theTjoundaries is given by the negative of expression (201) of 
reference 5. If to this part is added the potential energy acquired by the lateral loads, the resulting expression for the potential 
energy of the external forces is 

Equation (27) applies only when tlie reactions do no work and therefore acquire no potential encigy in the couree of the 
plate’s deflection. The most commonly assumed boundaries satisfying this condition are free, simply supported, and clampeal 
edges. The potential-energy expression for plates with more general boundary conditions must include terms correspomling 
to the work of the reaction forces. This more general case is considered in appendix C. 

In this section equation (27) has been established by means of physical reasoning. A more rigorous deiivation of equa- 
tion (27) for the special case of a rectangular plate is given in appendix D. 


POTENTIAL ENERGY OF SYSTEM 


The total potential energy V of the system comprising the plate and the forces acting on it is the sum of the strain energy 
Vi and the potentital energy of the external forces F* or, by addition of equations (26) and (27), 


F=iffi 

2j J h— DqJj 


Dxfiy+Dyiix r ^ /dw _ r ^ /'bw Qy X~j 

1— MxMy J^ax)j\j>y\^ ^/J 




1— u 


r a (m Qy \-l r d /dte Qy \ , b Qx* I Oil . 

ay U + 2 Lax w S7j+a^ 


(g) + i., (^)’+2^., ^ |] * 


(28) 


The above expression applies when the boundaiy reactions do no work and thei’efore acquire no potential energy in the 
course of the plate's deflection. This equation is therefore applicable when the edges of the plate are free, simply supported, 
or clamped. The potential-energy expression for a plate with more general boundary conditions is given in appendix C. 
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RECAPITULATION OF PRINCIPAL RESULTS 

1. The physical properties needed for smali-deflection 
analysis of an orthotropic plate in which deflections due to 
transverse shear are to be considered are the flexural stiff- 
nesses Dx and Dy, the corresponding Poisson ratios pj and 
/xy defined in terms of curvatures, the hvisting stiffness D^g, 
and the transverse shear stiffnesses Dq^ and Dq^. These 
constants can be evaluated theoretically or by. tests on 
samples of the plate as described in appendix A. Four of 
these constants are related by the reciprocal relationship 
HxDg=figDx derived in appendix B. 

2. The differential equations relating the deflections w, 
the lateral load q, and the internal forces and moments Nx, 

gf iV * fgf Qxj ^gj H-f*) ilTj,, and are 


3*10 Mi , 


bx 


bho_ Mg, Mi 1 bQg 

by*- Dg^f^^ Di-^DTby 


b*to Alig , 1 1 


. _ ^ 

bxby Big 2 by ~^2 Bq^ bx 


^Qt t 1 1 '''tc,g 
Xo, "1*0 71 X-, 


relating distortions to distorting moments and forces, and 

bMxg , bMt 


Qs 


Qg=- 

for equilibrium. 


by 

bflf„ 


bx 

bMg 


bx 




3. The first three equations can be solved for AIx, Mg, and 
Alxg to obtain 

M [~£> /c>io QA , 

^ l-iiiij.glpx\bx BaJ~^^''by\by 

If Dg ra /bw Qg\.^ 

* l—tixtigL^\^y BgJ'^^’^bx\bx T^Jj 

\f 1 71 r b /bto Qg \ ■ b /bw 

Substitution of these expressions into the last three equations 
and solution of the resulting equations by means of opera- 
tional determinants give the foUovring differential equations 
with variables separated, for the case in which Nx, Ng, and 
Nxg are constant throughout the plate: 

[B]w=-[M]q 

[B]Qx=-[N]q. 

[B\Qg=-[P]q 


where [B], [il/], [iV], and [F] are differential operators defined 
by equations (14). 

4. Three types of support commonly assumed at the 
boundaries of a plate are no support . (free edge), simple 
support, and clamping. These types of support can be 
described in teims of deflection, shears, and moments for an 
edge parallel to the y-axis as follows: 

For a free edge, 

Mx=Q 

Mxg=Q 

„ , buj , bio 

Qz+Nx £)a;+-^*y£)y — 0 

For a simply supported edge at which the support is applied 
over the entire thickness, 

10=0 
A/,=0 
Qg r. 

ll— — 0 


For a simply supported edge at which the support is applied 
only to the middle surface, 

10=0 

Mx=0 

Mxg=0 

For a clamped edge at which the support is applied over the 
entire thickness, 

10=0 

bio 

bx Ba^~^ 

-^—=0 

D,, “ 

For a clamped edge at which the support is applied only to 
the middle surface (a type of support very unlikely to be 
met in practice), 

10=0 

^ Q*._n 

bx B^^~^ 

Mxg=0 

The conditions for an edge parallel to the x-axis can be 
written by replacing x by y and vice versa, except in the 
subscripts of Alxg and Nxg. 

Boundary conditions can also be written for more general 
types of support. (See appendix C.) 
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5. The potential energy of a plate in which the middle- 
surface forces are assumed to remain unchanged in the 
course of the plate’s deflection and for which the moments 
and vertical forces at the boundaries do no work is 



The most important types of boundary to which this expres- 
sion applies are free, simply supported, or elamped. For 
more general types of support, in which the boundary reac- 
tions do work in the course of the plate’s deflection, the 
potential-energy expression must bo extended to include 
terms representing the potential energy' of the reactions. 

The calculus of variations can be used to show that iu 
order for the potential energy to be a minimiun the differ- 
ential equations of equihbrium and the boundary conditions 
must be. satisfied. (See appendix C.) 


Langley hlEMOEiAL Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., September SO, 1947. 


APPENDIX A 


TESTS TO DETERMINE PHYSICAL CONSTANTS 


The pui'pose of this appendix is to give descriptions of 
possible tests for determining the physical constants. No 
consideration has been given to details of testing technique. 
Practical considerations may dictate . changes in the test 
procediures described or the. quantities to be measured. 
These changes, however, willno tbeof f undamentalimportance. 

Test for I?, and — ^The flexural stiffness can be de- 

termined by cutting a beam from the plate in the a!-direction 
and loading it as sho%vn in figure 3. The supports and load- 
ing should be such as to make for minimum interference with 
the anticlastic curvature. The middle section is subjected 

only to a pine moment Pd; the cuiwature ^5 in this part 

can be determined from deflection or strain-gage measui’e- 
ments, and the flexural stiffness is given by equation (1) : 


L»x 




Pd . 
d^w 
' dP 



where b is the width .of the beam. If .the transverse curva- 

ture ^ is measured (the beam must be wide enough to 

permit accurate iheasurement of this emwature) , the Poisson 
ratio /ix can be calculated from equation (2) : 


d^ip 


-(A2) 



in figm’e 4. 'The beam at any station x is subjected to a known 
bending moment Mj equal to | ^ 

change of transverse shear equal to — The curvature 


d'w 

dr* 


along the beam can be determined from deflection meas- 


urements. (Sti'ain-gage measurements on up])cr and lower 
surfaces of beam are inappropriate because curvature due to 
rate of change of shear is not accompanied by stretching of 
the surfaces.) The flexural stiffness having been previ- 
ously determined and the transverse moment taken to 
be zero, equation (9a) can be solved for Dq^ to obtain 



Test for Dq ^. — The transverse shear stiffness Dq^ can be 
determined by loading the beam with a uniform load as shown 
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Fiquse 4.— Test to determine Da,. 


Tests for and — The constants I>„ ii^, and 

can be determined by tests similar to those already described 
but on a beam cut in the y-direction. 

Test for — ^The t^visting stiffness can be determined 

by cutting a rectangular panel from the plate, two edges 
parallel to the a;-axis and two edges parallel to the y-axis, 
placing some reinforcement at the edges to keep the 
boundary cross sections rectangular, and loading the 
panel at the comers as shown in figure 5. This loading is 
statically equivalent to a twisting moment AIx, distributed 

p 

around the edges and equal to -g' If the edge reinforce 

ments keep the boundary cross sections rectangular, then 
the shear angles y* and yg can be assumed to be zero and 



the plate to be in a condition of pure twnst with no 
or Qg loading present. The twist can be calculated 

from the measured comer deflections as 


yu? 45 
bxby iL 


The stiffness Z?,, is calculated from the fonnula that applies 
when only is acting, namely, equation (5) : 


Dxg — 


P/2 PIL 
'¥w~UllL 86 


i>xby 


(A4) 


APPENDIX B 

DERIVATION OF RELATIONSmP 


Betti’s reciprocal theorem (reference 10) can be expressed 
as follows: Let two groups of forces be applied to a struc- 
ture, each group of forcre producing distortions that me 
directly proportional to the magnitude of the forces; then, 
the work of the first group of forces acting through the dis- 
placements produced by the second group is equal to the 
work of the second group acting through the displacements 
produced by the first. 

The structure tn which this principle is applied is the ele- 
ment dx dy of figure 1. Let the first group of forces consist 
of the moments dy. The distortions produced are the 

cmwatures and where, from equations (1) and (2), 

ilii 
~ Dx 


and 


y-io Mx 


The second group of forces are the moments Aig dx, and the 
distortions produced by the group are the curvatures^- 

and ^ where, from equations (3) and (4), . 

W‘~ 


and 




AL 


Dg 


The work done by the first group of forces M* dy in associa- 

M 

tion with the curvature produced by the second 


group 15 


or 


-AIxAIg ^ dx dy 


Similarly, the work done by the second group of forces 

A£ 

Alg dx in association with the curvature /i* produced by 
the first group is 


—Alydx 


(y-z 


or 


—AlxAdg ^ dx dy 


Equating the expressions for the two W'orks and eliminating 
the common factor —AlxAig dx dy give 

l^g Dx 

from which is obtained equation (8). 



APPENDIX C 

DERIVATION OF EQUILIBRIUM EQUATIONS AND GENERAL BOUNDARY CONDITIONS BY A VARIATIONAL METHOD 

In the body of this paper only free, simply supported, and clamped edges were considered. These typos of ])oundary 
conditions are characterized by the condition that the moments and vertical forces at the boundaries acquire no potential 
energy as a result of the plate's deflection. This condition holds by virtue of the fact that either the moments and forces at the 
boundaries are zero or the points of application of the nonzero boundary reactions do not move. A more general type of sup- 
port, in which neither of tlicse conditions holds, is discussed in the following section. 

Potential-energy expression. — For simpheity a rectangular plate with edges a:=0,a and y=0,h is considered. The 
boundary reactions of the plate consist of distributed bending moments, twisting moments, and vertical forces statically 
equivalent to the limits of the internal moments and sheai^s boundaries are approached. The intensities of the reactions 
(moment or force per unit edge Icngtli) are denoted by Mz, Mxtn and "5* along those boundaries parallel to the y-axis and 
M*, ~Mzy, and Tly along those boundaries parallel to the z-axis. (Note that the symbols used for the reactions are distinguished 
from the corresponding symbols for the internal forces by means of bars placed above the symbols.) 

The potential energy of a plate the edges of which are other than free, simply supported, or clamped can be^ wTitten 
by adding to equation (28) line integrals representing the potential energy of the reactions. The resulting expression is 


y=i 
' 2 




AixAi»L3a;\2te IDqJj 


DxlJ-v+DslJ-z l 

l—tJ-zHy 


Ldr V&a- DqJ\ 1 —tixfiy ~ 


j.‘ [-». I : j: [-®- @-^>5-] I : (c*) 

In the last two integrals, representing the potential energy of the boundary shears and moments, and 

ai-e the rotations parallel to the zs-plane and ys-plane, respectively, of. an originally vertical line element in the edge of 

the plate.' 

Mininuzation of total potential energy. — The conditions that must be satisfied if the total potential energy V of the 
system is to be a minimum are now considered. By the calculus of variations (reference 11), minimization of T requires 
the vanishing of the first variation SV. The first variation can be evaluated from equation (Cl) as 



«; 2 Dz [~_d /^_ Qz\] 




dx 


+ 


\ l—fixlJty 




i>5 


Dxl^y+D,l^x \ r ^ Qx'NI 

. I— fixity Ju>y\^y ^9yJj 


d5 


2Dy r ^ /dw Y 

1— ^yJ- 


\(>X VqJ 


dx ^ dy 


d6 




dy 


2Q, 


./bv) 

Qy\ 


dJ 

dy 



^Qy) 

+ 

__ 



+ 


®^L^xVdy Z>Q^/'^dy\dx 

Those terms in the above expression that contain derivatives of integrated by 


parts so as to reduce the order of the derivatives. The resulting expr®sion for SV contains surface inh^rals of the type 

rr, Qx\^..,.. 
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(...)^dxdy 


and 


JJ 

— J* J* (...) S ^ dx dy 


Q. 


The first of these two integrals can be int^ated again by parts so as to leave only terms containing 5w rather than deriva- 
tives of 5w. With the aforementioned int^rations by parts performed and after considerable rearrangement of terms, the 
expression for SV becomes 


sir- r ^ D- QAi, 

Jo Jo h— MxM^L5^\2>a: ^Jy2\ 1— DqJj 

1 r_^ A M, D» r ^ 

2 \ 1 — J \j>a^ by \by Bq^ JJ'*" 1 — /ixM, W 


+ 


B. 


/boo_ 


1 1 

fbw 


_a~N ^~ 

-iV — - 

-2JV 1 

W~ 

“V 

' bxby*’ 

i^bx 

dJ. 



bxby] 


Jo Jo U— Ati;ivL2>2^ Bo J_y2\ l— M xM» /L2>a:b2/W DaJj 


i®*>[5l^^-^;)+$(S-^)]+e- ! * St^^y+ 

( 1 f B^.+D.y. \ r Qr \~|, D> - 

Jo Jo I 2 V l—Mxth J Lsia; ^ \c52 WjJt DqJS^ 


b* fbw 




B, rj!i J* fDdh±D^\ 


1 

('bU!_ 

<2 A 



i5o,/ 


l-y^tig\J)o^\bx SZJJ 2\ 1— /L^rdyVc>3/ 


- B 

o -^xv 


~ y /i 
_dxby \i 




[I ( 1 -^.)]+®'}* l>»+ 

( _ 1 / z>^,+z>,;,A r y /^_ \i_ x>, 

Jo ( 2\ 1— ;ixMr /L2>J:£)2/ \Sx 1— W £^o,/J 

1 n r2>’* S, \ , b* /bio QrM I Ar I AT- ^ 71 Is I* .7 i 

2^*'L&x"Cby I5^,)+bxb2/Cbx 


^ ( Bxliy-{-ByIl:\ f— ^• 


2V 1- 


MiMjr 




^ gxNl, I>, fd/bUT <2A1,V7 ) ’‘dxA- 
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By virtue of equations (12) and (8), the above expression for SF can bniewritten as .... 


5F= 


b^Mx 


b^My ( 

bz? 

^ bzby 

by‘ \ 


;o JO 


'-(s+iv, g+Af, ^+m, ^)] !«.&*+ 


:(■ 


j , . r* r* / aM, , 3M„ 

(--5^+ 


dx 


■Qy) 


sQv 


dz dy+ 






dx dy 


1-^-) * (s-4) !/■'+ 


j; a |>,+ j; g-5,) a» |>.- 

j; a 1;^.+ j; 

In order for 5V as given by the above expression to be zero for all possible values of Sw, 8 S^, and 5 the various 

^Qx ^Q\f 

integrals must individually be zero. The following differential equations result from equating the surface integrals to zcroi 


yii4 

„ b^Mxyj 

,^My. / 

bz^ 

bzby 

^ by^ \ 


^+2AT,, 


Q. 

Qv 


bMx^ ■ bMt 
' by ^ bx 

bMiu . bM^ 
bz by 

/'bMt bMt 


(C2) 


By virtue of the last two of equations (C2), ^ integrals can bo replaced by (?* 

and Qt, respectively. Equating the line integrals to zero then gives the following boundary conditions required to insure 
that 5F=0; 

At x=0,a 


At 3/==0,f 


n , TiT ^ , -KT ^ 

or 


- 

ew=o'| 



11 

or 

s ( 

f'btv_ 

^bx 

o 

11 

- - 

(C3) 

II 

or 

a 

/bw 

W 

\ 

O 

II 



^ ■ •%'r 1 Vt ^ 

by'^^^ bx~^'' 

or 

- 

- 

> 

aw=o 



My=My 

or 

a 

/bw 

O 

11 

1 

► 

(C4) 


or 

a 

/bw 

V&x 

A 

o 

II 

1 




Equations (C2) are the differential equations that must be satisfied if the potential energy is to bo a minimum. They 
wili be recognized as the equations of equilibrium, equations (11). 

Equations (C3) and ((Dd) are the boundary conditions that must be satisfied if the potential energy is to be a minimum. 
The left-hand groups of equations (C3) and (C4) imply that the limiting values of the internal forces and moments, as the edge 
of the plate is approached, must be in equilibrium with certain prescribed forces and moments externally applied at the edge 
(the prescribed forces and moments being designated by means of the horizontal bars). The right-hand groups of equations 
(C3) and (C4) imply that the displacements at the edge must have certain prescribed values. 

The boundary conditions given by equations (15) to( 19) for free, simply supported, and clamped edges parallel to the y- 
axis are special cases of equations (C3). For example, the boundary conditions for a simply supported edge (equations (16)) 

can be obtained from equation (C3) by prescribing the values of w, Mx, and boundaries z=0, a. 
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If a plate is elastically supported at the boundaries, the dastic support may sometimes be conveniently thought of as 
made up of three rows of dosdy spaced discrete springs at each edge: a row of deflectional springs, a row of rotational springs, 
and a row of torsional sp ring s, having the known stiffnesses per inch ki, ks, and ka, which may vary along the edge. For this 
type of support the vertical shear reaction at any point along the edge is proportional to the vertical deflection at that point 
and the twisting and bending moment reactions are proportional to the corresponding rotations of an originally vertical line 
element in the edge. The boundary conditions for this type of support can be obtained from equations (C3) and (C4) by set- 
ting at z=0 

Q.=hv, 

at x=a 
aty=0 

Q,=hw - kt 

-k,w Tl,=h 

The signs in the above boundary conditions foUow as a result of the directions assumed for positive shears and moments. 



APPENDIX D 


DERIVATION OF EQUATION (27) FOR THE POTENTIAL ENERGY OP THE EXTERNAL FORCES 


A rectangular plate the edges of which are x=0,a and 
y=0,b is considered (fig. 6). The boundary conditions as- 
sumed are the usual conditions corresponding to zero work 
by the reactions; that is, each edge is either free, simply 
supported, or clamped. 

The horizontal loads iV,, Nf, and are assumed first to 
be applied at the boundaries with no lateral load. As a 
result the middle plane (and aU horizontal planes) of the plate 
stretches; thus, the constant stretching energy discussed 
previously in coruiection with the strain energy of the plate 
is produced, and slight shifts in the points of application of 
the edge forces iV*, and JV*y are caused. These new 
positions of the points of application are used as the arbitrary 
fixed reference points in any future measurements of the 
potential energy of the horizontal edge forces. 

If the lateral load g is now applied, the middle surface 
acquires the displacements w(x,y) in the s-direction, u(x,y) 
in the ar-direction, and v(x,y) in the y-direction. As a result 
of these displacements, the lateral load acquires the potential 
energy 

~ Jo Jo 


and the edge forces acquire the potential energy 




+ 


NiyU 


;) 


dx 

(D2) 


The moments and vertical forces at the boundaries do no work 
and therefore acquire no potential energy during deflection. 



Fiovbe 6.— Eectangolar plate with hodzoatal forces applied to boundaries. 


By use of the formula for integration by parts, expression 
(D2y for the potential energy of the edge forces can be re- 
written in terms of the interior forces and displacements as 
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In the development of the differential equations and in this 
section the middle^surface stresses Nx, Ng, and Nx^ are as- 
sumed to remain .unchanged in the course of the plate’s 
deflection. Equations (10) for equilibrium of horizontal 
force, consequently, remain satisfied at all times, and, there- 
fore, the last two integrals of expression (D3) vanish. Further- 
more, the assumption that the middle-surface stresses remain 
unchanged implies that no stretching of the middle surface 
during deflection occurs. In order to prevent such stretch- 
ing the horizontal displacements u and v can be shown 
(p. 313, reference 6) to be related to the vertical displace- 
ments w as foUo\vs : . . 

ba: ™ 2 \ bz / 

^__1 
by~ 2\byJ 

bu . dv dw dw 
dy~^dx~ bz by — 

The fiist and only remaining integral of expression (D3) 
therefore becomes . 

3 XT S 1^] * 

Addition of expressions (Dl) and (D4) gives, as the total 
potential energy of the external forces, 

(D5) 

Although the derivation was carried out for the special 
case of a rectangular plate, equation (D5) also applies to a 


plate of any shape in which the middle-surface stresses 
remain unchanged during deflection. Equation (D5) is 
identical with equation (27). 
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